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Abstract. Using hydrodynamic simulations we investigate the rotational properties and angular momentum evolution of 
prestellar and protostellar cores formed from gravoturbulent fragmentation of interstellar gas clouds. We find the specific angu- 
lar momentum j of molecular cloud cores in the prestellar phase to be on average (j) = 7 x 10^" cm^ s ' in our models. This 
is comparable to the observed values. A fraction of those cores is gravitationally unstable and goes into collapse to build up 
protostars and protostellar systems, which then have (j) = 8 x 10" cm^ s"'. This is one order of magnitude lower than their 
parental cores and in agreement with observations of main-sequence binaries. The loss of specific angular momentum during 
collapse is mostly due to gravitational torques exerted by the ambient turbulent flow as well as by mutual protostellar encoun- 
ters in a dense cluster environment. Magnetic torques are not included in our models, these would lead to even larger angular 
momentum transport. 

The ratio of rotational to gravitational energy /? in cloud cores that go into gravitational collapse turns out to be similar to the 
observed values. We find, fj is roughly conserved during the main collapse phase. This leads to the correlation j oc M-'' , between 
specific angular momentum j and core mass M. Although the temporal evolution of the angular momentum of individual 
protostars or protostellar systems is complex and highly time variable, this correlation holds well in a statistical sense for a 
wide range of turbulent environmental parameters. In addition, high turbulent Mach numbers result in the formation of more 
numerous protostellar cores with, on average, lower mass. Therefore, models with larger Mach numbers result in cores with 
lower specific angular momentum. We find, however, no dependence on the spatial scale of the turbulence. Our models predict a 
close correlation between the angular momentum vectors of neighboring protostars during their initial accretion phase. Possible 
observational signatures are aligned disks and parallel outflows. The latter are indeed observed in some low-mass isolated Bok 
globules. 
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1. Introduction 

Angular momentum plays a pivotal role in star formation. The 
amount of specific angular momentum determines whether a 
collapsing protostellar core will form a single star or a binary 
or higher-order multiple system. 

Stars are thought to form by gravoturbulent fragmentation 
in interstellar clouds. The supersonic turbulence ubiquitously 
observed in molecular gas generates strong density fluctua- 
tions with gravity taking over in the densest and most massive 
regions. Once gas clumps become gravitationally unstable, col- 
lapse sets in and the central density increases until a protostellar 
object forms and grows in mass via accretion from the infalling 
envelope. Various aspects of the relation between supersonic 
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This dynamic picture of gravoturbulent star formation chal- 
lenges the so called "standard theory" where stars build up 
from the "inside-out" collapse of singular isothermal spheres, 
which are generally assumed to result from the quasistatic con- 
traction of magnetically supported cloud cor es due to ambipo- 
lar diff^usion tShu 1977: Shu. Adams. & Lizano 1987) . This 
picture, however, has always received strong criticism (e.g., 
'Whit worth & Su mmers" 1985': 'Whitworth et al."l996': 'Nakanol 
.19981, for a critical discussion). In particular, it seems strongly 
biased toward the formation of single stars dWhitworth et alJ 
1 19961) which is in contradiction to the observational fact 
that most (if not all) stars form as members of binary 
or higher-order multiple systems (see, e.g., the reviews by 
'Bode nheimer et al.ll2000t iMathieu et alJl200d and references 
therein). 
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Gravitational collapse in the astrophysical context always 
involves solving the angular momentum problem. It results 
from the blatant discrepancy between the specific angular 
momentum observed in low-density gas on large scales and 
the amount of rota tion present after collapse (Spi tzeJfl968t 
lBodenheimeJl995 ). The source of angular momentum on large 
scales lies in the differential rotation of the galactic disk and, 
closely related to that, on intermediate to small scales it re- 
sults from the high degree of vorticity inextricably adherent 
to turbulent flows. The typical specific angular momentum j 
of molecular cloud material, e.g., on scales of about 1 pc is 
cm s , while on scales of cloud cores, say below 
0.1 pc, it is of order of 10^' cm^ s A binary G star with a 
orbital period of 3 days has j x lO'^cm^s"', while the spin 
of a typical T Tauri star is a few xlO'^cm^ s"'. Our own Sun 
rotates only with j x 10'^ cm^ s"' . That means, during the pro- 
cess of star formation most of the initial angular momentum is 
removed from the collapsed object. 

The presence of magnetic fields, in principle, provides a 
viable mechanism for locally reducing the angular momentum 
through magnetic braking. This was treated approximately by 
lEbert. von Hoemer. & Temesvarv ( 1960), a nd lat e r calc ulated 
accuratelv bv iMouschovias & Pa leologoui ( 1 19791 Il980l) . The 
criterion for effective braking is essentially that the outgoing 
helical Alfven waves from the rotating cloud have to couple to 
the ambient medium over a volume that contains roughly the 
same mass as the cloud itself. For the strong magnetic fields 
required by the standard theory of star formation, the decelera- 
tion time can be less than the free-fall time, leading to efficient 
transfer of angular momentum away from collapsing cores, 
and thus, to the formation of single stars. Field strengths small 
enough to allow for binary formation cannot provide support 
against collapse, thus pointing toward a more dynamic picture 
of star formation as offered by gravoturbulent fragmentation. 

It is therefore a crucial test for any theory of star for- 
mation whether it can produce the required angular momen- 
tum loss during collapse while at the same time explain the 
high numbers of binaries and multiple stellar systems obs erved 
(e.g.. lDuauennov & Mavor"l99l': Hal bwachs e t al. '20031 ). In a 
semiempirical analysis of isolated binary star formation Fished 
presented the effects of turbulence in the initial state 
of the gas on binary orbital parameters. These properties were 
in agreement with observations if a significant loss of angular 
momentum was assumed. In the current investigation we fo- 
cus on models of non-magnetic, supersonic ally turbulent, self- 
gravitating clouds and analyze the time evolution of angular 
momentum during formation and subsequent collapse of pro- 
tostellar cores. Our main question is whether gravoturbulent 
fragmentation can solve or at least ease the so called "angular 
momentum" problem without invoking the presence of mag- 
netic fields. 

The structure of the paper is as follows. In Section|2lwe in- 
troduce and discuss the numerical method to calculate the dy- 
namical cloud evolution and the suite of models included in the 
current analysis. In Section |3l we present results on the angu- 
lar momentum distribution of starless molecular cloud cores. 
We call them prestellar cores. Some of them collapse to be- 
come protostellar cores. In Section |3 we investigate the angu- 



lar momentum evolution of their central protostellar objects. 
We report a statistical correlation between specific angular mo- 
mentum and mass, and analyze its dependence on the turbu- 
lent environment. The angular momentum vectors of neighbor- 
ing protostars tend to be aligned, at least in the early accretion 
phase. This is discussed in Section|5l Finally, in Section|6l we 
summarize and conclude. 



2. Models of dynamical cloud evolution 

2.1. Numerical method 

To adequately describe the angular momentum evolution dur- 
ing gravoturbulent fragmentation, it is prerequisite to resolve 
formation and subsequent evolution of collapsing fragments 
over several orders of magnitude in density. Due to the stochas- 
tic nature of supersonic turbulence, it is not known in advance 
where and when local collapse occurs. We therefore resort to 
smoothed particle hydrodynamics (SPH) to solve the equations 
of hydrodynamics. It is a Lagrangian method with the fluid rep- 
resented by an ensemble of particles and flow quantities ob- 
tained by local averaging (.Ben z 1990; Monag han 1992) . The 
method is able to resolve large density contrasts as particles 
are free to move and so naturally the particle concentration 
increases in high-density regions. Because it is computation- 
ally prohibitive to treat the cloud as a whole, we concentrate 
on subregions within the cloud and adopt periodic boundary 
conditions (Klessen 1997). Once the central region of a col- 
lapsing gas clump excee ds a density contrast of ~ 10^ , we in- 
troduce a "sink particle" teate, Bonnell, & Price" 1995*), which 
accretes gas from its surrounding while keeping track of mass 
and linear and angular momentum. Replacing the high-density 
cores with "sink particles" allows us to follow the angular mo- 
mentum evolution of collapsed cores over many free-fall times. 
Altogether, the performance and convergence properties of the 
method are well understood and tested against analytic models 
and other numerical schemes in the context of turbulent super- 
sonic astroph ysical flows. For a detai l ed discussion th e reade r 
is referred to 'Bate & Burkert' !l99D, 'Mac Low et al.' h99^, 
Lombai-di et il99 i), Klessen & Burkert (.2000, ,2001.) and 
Klessen et all ( I200(]Ik 



2.2. Models of turbulent self-gravitating clouds 

The large observed linewidths in molecular clouds imply the 
presence of supersonic velocity fields that carry enough en- 
ergy to count erbalance gravity on global sca les (see, e.g., 
the review by Williams. Blitz. & McKee 2000). However, it 
is known that turbule nt energy dissipates rapidly, roughly on 
flie free-fall timescale JMac Low et alJl998tlStone et alJl998t 
IPadoan & N ordlun? 1999). Unlike previously thought, this is 
independent of the presence of magnetic fields. Magnetic fields 
also do not significa ntly alter the efficiency of local coll apse 
for driven turbulence ( iHeitsch. Mac Low. & Kles sen'2001*). To 
a first approximation, the fields are therefore taken as being 
dynamically unimportant in our current models, and are not 
included. We also disregard possible feedback effects from 
the star formation process itself (like bipolar outflows, stellar 
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winds, or ionizing radiation from new-born O or B stars). Our 
analysis of angular momentum evolution of protostellar cores 
therefore focuses on the process of gravoturbulent fragmenta- 
tion, i.e. on the interplay between turbulence and self-gravity 
only. 

The suite of models considered here consists of 12 numer- 
ical simulations where turbulence is maintained with constant 
root-mean-square Mach numbers in the range 2 < Al < 10. 
This roughly covers the range observed in typical Galactic 
molecular clouds. We apply a non-local scheme that inserts 
energy in a lim ited range of wavenumbers at a given rate 
llMacLowlll999l) . We distinguish between turbulence that car- 
ries its energy mostly on large scales, at wavenumbers 1 < 
^ < 2, on intermediate scales, i.e. 3 < A: < 4, and on small 
scales with 7 < ^ < 8. The corresponding wavelengths are 
{ - L/k, where L is the total size of the computed volume. The 
models are labeled mnemonically as MAikk, with rms Mach 
number Ai and wavenumber k. Each of these simulations con- 
tains 205379 SPH particles. We also consider a model that 
is globally unstable and contracts from G aussian initial con- 
ditions without turbulence (for details see iKlessen & Burkert' 
boOO, 2Q0D. It is called GA and was run with 500000 parti- 
cles. The main parameters are summarized in Table [J Note, 
the final star formation efficiency varies between the different 
models, as indicated in Column 5 of Tabled This simply re- 
flects the evolutionary stage at the time when we stop the cal- 
culation. In some cases the accretion timescale is too long to 
follow the simulation to high efficiencies. 

Table 1. Sample parameters, name of the environment used in 
the text consisting of the Mach number M and the driving scale 
k (GA denotes the model with Gaussian density), number N of 
protostellar objects (i.e. "sink particles" in the centers of proto- 
stellar cores) at final stage of the simulation, percentage of ac- 
creted mass at final stage Macc/M,o,, parameter ^ see Eq. (|4}, 
parameter S see Eq. (|8} 



Name 


k 


M 








S 










[%] 


flO^o 


cm^ s^'j 


M2.0k2 


1..2 


2.0 


68 


75 


1.7 


1.7 


M2.0k4 


3..4 


2.0 


62 


48 


2.0 


2.0 


M2.0k8 


7..8 


2.0 


11 


66 


1.7 


1.6 


M3.2k2 


1..2 


3.2 


62 


80 


1.2 


1.3 


M3.2k4 


3..4 


3.2 


37 


82 


2.0 


1.9 


M3.2k8 


7..8 


3.2 


17 


60 


2.7 


2.6 


M6k2 


1..2 


6.0 


110 


76 


1.3 


1.3 


M6k4 


3..4 


6.0 


60 


65 


1.5 


1.7 


M6k8 


7..8 


6.0 


7 


4 


1.9 


1.4 


M10k2 


1..2 


10.0 


100 


38 


1.0 


1.0 


M10k4 


3..4 


10.0 


10 


6 


2.0 


1.4 


MlOkS 


7..8 


10.0 


27 


8 


1.4 


1.05 


GA 






56 


85 


1.4 


1.05 



2.3. Physical scaling and naming convention 

The models are computed in normalized units using an isother- 
mal equation of state. Scaled to physical units we adopt a 



temperature of 11.3 K corresponding to a sound speed Cj = 
0.2 kms"', and we use a mean density of n(H2) = 10^ cm"^, 
which is typical for s tar-forming molecu lar cloud regions (e.g. 
in p-Ophiuchus, see iMotte. Andre. & NerL.1998J . The mean 
thermal Jeans mass' in all models is (Mj) = 1 Mq. The tur- 
bulent models contain a mass of 120 within a cube of size 
0.28 pc, and the Gaussian model has 220 Mq in a volume of 
(0.34 pc)-*. The global free-fall timescale is = 10^ yr, and 
the simulations cover a density range from «(H2) ~ 100 cm"-* 
in the lowest density regions to n(H2) ~ 10^ cm"^ where the 
central parts of collapsing gas clumps are converted into "sink 
particles". 

Each "sink particle" defines a control volume with a fixed 
radius of 560 AU. We cannot resolve the subsequent evolution 
in its interior After ~ lO^* yr a protostar will form in the very 
center Because of angular momentum conservation most of 
the matter that falls in will assemble in a protostellar disk. 
It is then transported inward by viscous and possibly grav- 
itational to rques (e.g., 'Bode nheimeJ 1 1 9951 IPapaloizou &Linl 
ll995tlLiir & Papaloizou 1996). With typical disk sizes of order 
of several hundred AU, the control volume fully encloses both, 
star and disk. If low angular momentum material is accreted, 
the disk is stable and most of the material ends up in the central 
star In this case, the disk simply acts as a buffer and smooths 
eventual accretion spikes. It will not delay or prevent the mass 
growth of the central star by much. However, if material that 
falls into the control volume canies large specific angular mo- 
mentum, then the mass load onto the disk is likely to be faster 
than inward transport. The disk grows large and may become 
gravitationally unstable and fragment. This will lead to the for- 
mation of a binary or higher-order multiple ( iBodenheimer et all 
2000). 

Throughout this paper, we adopt the following naming con- 
vention: In the pre-collapse phase, we call high-density gas 
clumps prestellar cores or simply gas clumps. They build up at 
the stagnation points of convergent flows. The flows result from 
turbulent motion that establishes a complex network of inter- 
acting shocks. The gas clumps are identified and characterized 
using a three-dim ensional clump-find algorithm as described 
in Appendix A of ' Klessen & BurkertI llOOd) . The fluctuations 
in turbulent velocity fields are highly transient. They can dis- 
perse again once the converging flow fades away. Even clumps 
that are strongly dominated by gravity may get disrupted by 
the passage of a new shock front. For local collapse to result in 
the formation of stars, Jeans-unstable, shock-generated density 
fluctuations therefore must collapse to sufficiently high densi- 
ties on time scales shorter than the typical time interval be- 
tween two successive shock passages. We include in our analy- 
sis only Jeans-unstable gas clumps. Angular momentum is cal- 
culated from the internal motions with respect to the location 
of the density maximum. These objec ts correspond to the so 
called starle ss cores observed, e.g., by Goodma n et al 
iBarranco & Goodman ( 1998), Jijina et al. ( 1999), and others. 

' We use a spherical definition of the Jeans mass, Mj = 4/3 npA^, 

with density p and Jeans length A; = (^^) ' and where G and K are 
the gravitational and the gas constant. The mean Jeans mass (Mj) is 
then determined from the average density in the system (p). 
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They are thought to collapse and build up a central protostar or 
protostellar system in the later stages of evolution. Once col- 
lapse has lead to the formation of an embedded protostar (in 
our scheme identified by a central "sink particle") we call the 
object protostellar core or also protostar. The angular momen- 
tum is obtained as the spin accumulated by the "sink particle" 
during its accretion history. The distribution is best compared 
with observations of main-sequence binaries as we expect the 
unresolved star-disk system interior to the "sink particle" to 
break up into a binary or higher-order multiple. 

3. Molecular cloud clumps and prestellar cores 

Figure ^ shows the distribution of the specific angular mo- 
menta of the gas-clumps that were identified in the turbu- 
lent environment M6k2 (see Table We compare with 
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Fig. 1. The distribution of specific angular momenta of prestel- 
lar cores formed in our simulations using model M6k2 (non- 
hatched histogram) is compared to the distribution of spe- 
cific angular momenta of observed molecular cloud cores 
(hatched distributio ns). The obse r vation al data were taken in 
(a) from Table 5 i n Caselli et al. (2002), in (b) from Table 2 
in iGoodman et alj ( .1993.). Table 4 in Barranco & GoodmanI 
(Il998l) an d Table A2 in Ijiiina et all t\99^ and in (c) from 
Table 7 in IPirogov et al.l ( l2003l) . ^fe take fj to represent the 
percentage of the total number of existing cores in a certain 
specific angular momentum bin. 



observational values taken from Caselli et al. ' hmf) for (a), 
from iGoodman et alJ (il993.) . .Barranco & Goodmaa C1998,) . 
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Fig. 2. Distribution of specific angular momenta of the proto- 
stars or protostellar systems at different evolutionary phases of 
the numerical model M6k2 as denoted by the local star for- 
mation efficiency in percent ((a)-(d)). We compare in (e) with 
the /'-distributi on of binaries among nearby G-dwarf stars from 
IpuQuennov & MavoJ ( Il99ll) (for details see text) and in (f) 
with the distribution of specific angular momenta of binaries in 
the Taurus star-forming region from Figure 5 in .Simoa 
Again, fj represents the normalized distribution function. 



Ijiiina et al.l(ll999l) for Tb') and fromlPiroeov et alJ ( Eooi for (c). 
According to Goodman et al. '("1993^ the values for the specific 
angular momenta are obtained by using best-fit velocity gra- 
dients from maps of observed line-center vel ocities under the 
assumptio n of solid b ody rotation. The cores in'Goodman e t alJ 
Il993). Barranco & Goodman ( 1998) . Jiiina et al. ( 1999) were 
mapped in the ( J.K) = ( 1 , 1 ) transition of NH ^, whereas the 
massive cloud cores in Pirogov et al. (2003) and the dense 
cloud cores in Caselli ^ et al. ( 2002 ) were mapped in N2H^. 
As found in Casell i et al.l (l2002l)~ the two molecular species 
trace essentially the same material, especially in starless 
cores. Model clumps and observ ed molecular cl oud clumps 
have com parable mass spectta dKlesseij |2001ai) and simi- 
lar shapes essen & Burkert'2000':'Ballesteros-Pare des et alJ 
12003). Fi gure [2 demonstrates that similar holds for the distri- 
bution of specific angular momenta. 

Note, when transforming from dimensionless code units 
into physical scales, the specific angular momentum depends 
on the mean density n and the temperature T as j cc T/ Vn(H2). 
In Figure QJ, we use our standard scaling corresponding to re- 
gions like the p-Ophiuchi main cloud tMotte et al..l998a) . This 
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is ade quate for the low-mass cores studied by ICaselli et alJ 
1 I2OO2I) . and in Sect. |3 we will furthermore show that the spe- 
cific angular momenta of collapsed cores then fall into the right 
range for main-sequence binaries. 

We find that the specific angular momentum of prestellar 
cores has values between 1 x 10^" cm^ s"' and 5 x 10^' cm^ s"' 
with a mean value of approximat ely (/') ^ 5 x 10^° cm^s~'. 
This is in good agreement with the lCaselh et al.l ( l2002h sample 
which has {j} = 7 x 10^"cm^s"'. Their cloud cores have a 
mean mass of ~ 6 Mq comparable to the core masses in our 
simulations. A Kolmogorov-Smirnov (KS) test was performed. 
We find that at a 50 % level the distributions are statistically 
indistinguishable. 

However t he cor es observed bv lGoodman et alJ ( Il993l) and 
IPiroeov et al.l (l2003h trace lower densities and have higher 
mean masses of around 50 Mg and 500 Mq, respectively. In 
Figure^ we therefore use n(H2) = 1 x 10"^ cm"^ and T = lOK 
leading to (/' ) - 1.5 x 10^' cm^ s~'. This matches the the obser- 
vations since lGoodman et alJlll993|) find (j) - 6x 10^' cm^ s"' 

and a medi an value of 3 xlO^^cm^s"'. The massive cores 

• 1 I 1 

mapped by Pirogov et al. (2003) have higher velocity disper- 
sions, higher kinetic temperatures (20 - 50 K) and a den- 
sity n(H2) > 1 X 10"* cm""*. The resulting mean specific angular 
momentum is {j} = 1.5 x lO^^cm^s"'. Again, adequate scal- 
ing in Figure^ results in higher values of j in our models and 
leads to better agreement with the observed distribution. 

Given the simplified assumptions in our numerical mod- 
els, we altogether find remarkably good agreement with the 
observed specific angular mo menta in the prestella r phase. 
Similar findings are reported bv Gam mie et alJ ('2003). Similar 
to our study, they follow the dynamical evolution of isothermal, 
self-gravitating, compressible, turbulent ideal gas. However, 
they include the effects of magnetic fields and solve the equa- 
tions of motion using a grid-based method (the well-known 
ZEUS code). Their approach is thus quite complementary to 
ours. The /' distribution that results from their simulations 
peaks at 4 X lO^^cm^s"'. They fix the mean number density 
at n(H2) ~ 1.0 X 10^ cm^^ and use T = 10 K. Using the 
same physical scaling we get very similar values, i.e. {j} = 
4x lO^^cm^s-i. 

This mean value also falls into the range of spe cific an- 
gular momenta of cores that form in simulations by Li et al.' 
12004). They also use a version of the ZEUS code (ZEUS-MP) 
to perform high-resolution, three-dimensional, super- Alfvenic 
turbulent simulations in order to investigate the role of mag- 
netic fields in self-gravitating core formation within turbulent 
molecular clouds. A dopting the same physical scaling as in 
iGam mie et al.' ('2003*), the specific angular momentum of their 
cores takes values between 5x10^' cm^ s"' and SxlO^^cm^ s 

4. Protostars and protostellar systems 

Figure 121 shows the distribution of specific angular momentum 
of collapsed cores at four different stages of mass accretion, 
ranging from 15 % of the total available mass in the molecu- 
lar cloud accreted onto collapsed cores in Figure 12^ to 60 % 
in Figure IJJl. While the distribution narrows during the evolu- 
tion, the mean specific angular momentum remains essentially 



at the same value 7 = (8 + 2) x 10'^ cm^ s"' with a range from 
10'^ cm^ s"' to 5 X 10^° cm^ s"'. The specific angular momen- 
tum of the protostellar cores in the considered model are ap- 
proximately one order of magnitude smaller than the ones of 
the Jeans-unstable clumps, but both distributions join without 
gap. In a statistical sense, there is a continuous transition as loss 
of angular momentum occurs during contraction. The range of 
specific angular momenta of the protoste llar cores agree wel l 
with the observed values for binaries ('e.g..'B odenheimeJl995l) . 
For this reason we compare in Figures |2t-d the model distri- 
butions with observat i ons of binaries among G-dwarf stars by 
IPuQuennov & MavoJ ( Il99ll) in Figure |2^ and with observa- 
tio ns of young s tar binaries in the Taurus star forming region 
bv lSimorJd 19921) in F igurel2f. 

IPuquennov & Mavoj (1199 li) derived a Gaussian-type pe- 
riod distribution for their sample. Based on this distribution we 
calculated the distribution of the sp ecific angular momenta us- 
ing the following equation (see also lKroup all995 b. Eq. 10): 

j = 6.23 X 10>« (1 - e^f^' pU^^HHai^ (em^ ,-^) 
^ ' (mi + mif^^ ^ ' 

where masses are in Mq and P in days. We used a primary mass 
m\ of 1 M0, a mean mass ratio between primary and secondary 
q - m^jmx - 0.25 and a mean eccentricity e - 0.31. 

The resulting distribution has a mean specific angular mo- 
mentum of 1 .6x10^*^ cm^ s"' , which agrees well with the values 
from the simulations. T his also holds for Figure|2f which was 
taken from Figure 5 in ISimonI (Il992h . The mean specific an- 
gular momentum here has a value of 1.6 x 10^"cm^s"'. The 
agreement of the distributions was confirmed by ■s.x^ statistical 
test. Since our numerical resolution is not sufficient to follow 
the subfragmentation of collapsing cores into binary or higher- 
order multiple systems, the time evolution of j is an important 
tool to evaluate or our models. We see a clear progression from 
the rotational properties of gas clumps (as discussed in Sect.|3} 
to those of the resulting collapsed c ores. A similar corr elation is 
observed between the cloud core (Caselli et al. 2002) and typ- 
ical main-sequence binaries ( Duauennov & Mavor.l991i) . The 
former may be the direct progenitors of the latter. 

4A. Example of the angular momentum evolution of a 
protostellar core in a cluster 

The evolution of the specific angular momentum of individual 
protostellar cores can be very complex depending on the rota- 
tional properties of their environment. There is a strong con- 
nection to the time evolution of the mass accretion rate. 

In Figure |3l we select five collapsed cores in model M6k2 
with about the same final mass. All of them show a similar 
evolution the specific angular momentum with increasing mass 
(Fig-Et) and time (Fig.|3j'). Nevertheless there are visible dif- 
ferences in the details. 

Initially, the specific angular momentum increases with 
growing mass. However, at later stages the evolution strongly 
depends on the secular properties of the surrounding flow. 
In cores 43, 96, and 101, for example, j decreases again af- 
ter reaching a peak value, while for cores 17 and 47 j stays 
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close to the maximum value. Depending on the specific an- 
gular momentum of the accreted material the resulting proto- 
stellar disks are expected to evolve quite differently. For ex- 
ample, preliminary 2-dimensional hydrodynamic calculations, 
show that core 17 and core 47 will probably develop a stable 
disk (Bodenheimer 2003, priv. comm.). The ratio of rotational 
to gravitational energy for the peak value of j is /3 - 0.005 
for core 17 and /3 - 0.003 for core 43. On the other hand, 
core 101 will fragment into a binary star. It has a y6 = 0.016. 
Also, the disk of core 96 is highly unstable with corresponding 
j6 = 0.016. The evolution of core 43 has not yet been followed 
sufficiently long to determine whether it will fragment to form 
a binary star or not. These results show the importance of the 
specific angular momentum on the evolution of the protostellar 
object. A high value /3 > 0.01 leads to fragmentation whereas 
B < .01 results in a stable disk. A similar result was found by 
IBossI (IT999) for slowly rotating, magnetic clouds. 

Figure (SJi shows, that the change in specific angular mo- 
mentum is closely linked to the mass accretion. At the point 
in time where the mass accretion rate (dotted line) has a pro- 
nounced peak, the change in specific angular momentum is 
also significant. A high mass accretion rate can result in an in- 
crease of specific angular momentum. But as seen for core 96 
in FigurelSj) a high mass accretion rate can also lead to a reduc- 
tion of specific angular momentum. The exact evolution of the 
specific angular momentum is thus closely linked to the flow 
properties of the surrounding material. 

4.2. Statistical correlation between specific angular 
momentum and mass 

As mentioned above we find a correlation between mass M and 
specific angular momentum j in a statistical sense. The result is 
depicted in Figure|3 It shows the angular momentum evolution 
as function of mass for all 1 10 collapsed cores in m odel M6k2. 

In Figure |^, following Goodman et alJ lll993h . we adopt 
rigid body rotation with constant angular velocity Q and and 
uniform core density p. With these assumptions the specific an- 
gular momentum j can be written as: 



(1) 



For a uniform density sphere = f . The mass M of a sphere 
with constant density po is related to the radius R via: 



M = —poR' 

From Eq.^and Eq.|2]follows that j can be expressed as: 



Anpo 



2/3 



2/3 



(2) 



(3) 



Therefore we fit the average angular momentum with a func- 
tion of the form: 



= J?[(M/Me) 



2/3 



(4) 



(3M \^/-^ rn 
. From Figure |4p the constant has 

a value of (1.3 + 0.6) x lO^^cm^s This fit formula can be 




Fig. 4. Absolute values of specific angular momenta {dotted 
lines) of all protostars (i.e. "sink particles") from our model 
M6k2 as a function of mass. The specific angular momenta 
are averaged at certain mass values which are separated by 
0.01 Mq and the resulting points are indicated by crosses. 
The solid line represents a fit of these averaged specific an- 
gular momenta in the mass range between and 1 .7 Mq In 
(a) we fit with a function of the form: ; = ^(M/Mq)^^^, 
= (1.3 ± 0.6) X 10^*'cm-s"', and in (b) we use a square 
root function: j = S y/M/Me, S = (1.3 ± 0.6) x lO^^cm^s"'. 
One standard deviation is marked by the dashed lines. 



applied to different turbulent cloud environments, and we list 
the corresponding values of ^ for our model suite in Tabled 
Using the fitted value and the density po = 4 x 10"'^ gcm"^ 



where protostellar cores are identified, we calculate an angular 
velocity Q = 1.33 x 10"" s"'. 

In this picture the ratio of rotational to gravitational en- 
ergy j6 can be written as: 



/3- 



(1/2)/Q2 
qGM^/R 



3p 

Snq pqG ' 



(5) 



where the moment of inertia is given by / = pMR^, and <7 = ^ is 
defined such that qGM^/R represents the gravitational potential 
of a uniform density sphere. With the assumptions of constant 
angular velocity and uniform density it follows that /5 is also 
co nstant. With Q and po as above, we get values 13 ^ 0.05. 

iGoodman et alJ ( Il993l) as well as iBurkert & BodenheimeJ 
teooo^ derived scaling relations where (3 is independ ent of ra- 
dius. Similar values for yS were also found by .Goodman et alJ 
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Fig. 5. Distribution of p obtained from our model M6k2 with 15 % (a), 30 % (b), 45 % (c) and 60 % (d) of the material accreted. 
For the collapsed protostellar cores in our models {solid lines) we calculate jS using Eq.|5j where we assume solid-body rotation 
of a uniform density sphere and take the calculated specific angular momentum j. For the prestellar cores (i.e. Jeans-unstable 
gas clumps; see the dashed lines) we derive p self-consistently from the three-dim ensional density and ve locity structure (using 
Eq. |9jl. For comparison we also indicate with dotted lines the values reported bv lGoodman et al] \ \99'3 . see their Figure 11) 
for observed molecular cloud cores. These were obtained with the same assumptions, and we use the same binning with fp 
representing the fraction of objects per p bin. 



for the observed cloud cores. In good agreement with 
our calculations they found that all values are below 0.18 with 
the majority under 0.05. 



The fit in Figure]^ rests on the assumption of the collapse 
of an initially homogeneous sphere with constant angular ve- 
locity. Using "sink particles" however, which have a constant 
radius, makes it necessary to examine another possibility. In 
Figure|4j3 we thus assume a constant radius R and a constant p. 
Choosing a constant p is supported by the observations as dis- 
cussed above and by our simulations as we show below. With p 



and Eq.|5l(which still holds) it follows that the angular velocity 
depends on the density as 



Q. = 



87rg 



Gpp. 



(6) 



Thus, the angular velocity D. is no longer a constant. This im- 
plies for the specific angular momentum j: 



(7) 



Following Eq.0we fit our data in Figure|4]3 with a square root 
function: 



J 



.!B^JM/Ms 



(8) 
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where the moment of inertia is given by / = pMR?, and q - j 
is defined such that qGM^/R represents the gravitational po- 
tential of a uniform density sphere. With the above density and 
angular velocity Eq.|5lresults in j3 ^ 0.05. For the new fit we 
find a scaling factor S = (1 .3 + 0.6) X 10^" cm^ s"' in the mass 
range < M/Mq < 1.7. Again, this fit formula can be applied 
to different turbulent cloud environments, and we list the cor- 
responding values of S for our model suite in Tabled 

The question remains if our simulations support the as- 
sumption of a constant /3. In Figure |5] we compare the distri- 
bution of /3 measured by Goodman et al. ( 1993) with values we 
extract from our model. For the prestellar cores, we use the 
definition 

(9) 

-'^grav 

and calculate rotational and potential energy, E^oi and figrav, 
consistently from the full three-dimensional gas distribution of 
each gas clump. We do not adopt any assumption about sym- 
metry and shape of the density and velocity structure. If we 
assume each clump is spherical and has roughly constant den- 
sity, as implied for example by Eq. |5] then /3 is overestimated 
by a factor of 2.7. This shows the importance of taking the 
full three-dimensional clump structure into account when ana- 
lyzing the rotational properties of molecular cloud cores. Both 
prestellar cores (i.e Jeans-unstable gas clumps) as well as pro- 
tostellar cores in our model typically have /3 < 0.3 with similar 
distributions. Thus, /3 stays mainly in the interval [0, 0.3] and 
in this sense it remains approximately constant during the col- 
lapse. 

It should be noted in passing, that we also looked at the den- 
sity structure of purely hydrodynamic turbulence, i.e. without 
self-gravity. If we again perform a clump-decomposition of the 
density structure and compute hypothetical y6-values, we find 
/3 ~ 2. This is indicative of the high degree of vorticity inher- 
ent to all turbulent flows. However, it also suggests that dense 
molecular cloud cores are strongly influenced by self-gravity. 
The fact that all cores in the observational sample have /3 < 0.2 
implies that gravitational contraction is needed to achieve den- 
sity contrasts high enough for sufficiently low /3 . This agrees 
with the picture of gravoturbulent fragmentation where molec- 
ular cloud structure as whole is dominated by supersonic tur- 
bulence but stars can only form in those regions where gravity 
overwhelms all other forms of support. 

Comparing the two fits in Figure |3 shows that our first set 
of assumptions is a better representation of the data. This is es- 
pecially true during the early accretion phase where we have 
good statistics. And it applies to different turbulent cloud en- 
vironments, as well. We conclude that - in a statistical sense - 
the angular momentum evolution of collapsing cloud cores can 
be approximately described as contraction of initially constant- 
density spheres undergoing rigid body rotation with constant 
angular velocity. This is consistent with the fact that cores from 
gravoturbulent fragmentation follow a Bonnor-Ebert-type ra- 
dial density profile ( Ballesteros-Paredes et al. 2003) and have 
roughly constant density in their in nermost regions , it also sup- 
ports the assumptions adopted by iGoodman et alJ ll 19931) and 
burkert & Bodenheimer.(.200a) . 



4.3. Dependence of the specific anguiar momentum 
on tlie environment 

When we compare the results of our complete suite of numeri- 
cal models (see Tabled we find as a general trend that the av- 
erage angular momentum falls with increasing Mach number. 
This is illustrated in Figure [6^. However, this follows mainly 
from a positive dependence of angular momentum on mass and 
from a correlation between average mass of the cloud cores and 
Mach number. As seen in Sect. l4.2l the specific angular momen- 
tum increases on average as the mass of the core rises. In envi- 
ronments with a low Mach number the mass growth of the cores 
is undisturbed over longer periods of time and so larger masses 
can accumulate. This can be inferred from Tabled where we 
list both the number of cores and the accreted mass. For higher 
Mach numbers more cores with on average less mass form. 
Thus, the average angular momentum is expected to decrease 
with increasing Mach number (Fig.|6t). 

To detect a direct dependence of the specific angular mo- 
mentum on the Mach number we select cores that belong in a 
certain mass bin and average the specific angular momentum 
only over those cores. The results are shown in Figures |6j3-|6ll. 
We find that there is in general little spread of specific angular 
momentum for different scales and different times in accretion 
history independently of the Mach number Low mass proto- 
stars (Fig.l^J)) are an exception, in low Mach number environ- 
ments they show especially low angular momentum in the early 
accretion phase. However, within the error bars we do not find 
a further dependence of j on the Mach number. 

Compared to the Gaussian collapse case, turbulent driving 
with small Mach numbers results in higher specific angular mo- 
menta. This is due to input of turbulent energy that can be con- 
verted into rotational energy if the turbulent velocities are not 
too high. 

5. Orientation of angular momentum vector 

We find in our simulations that neighboring protostellar cores 
have similarly oriented angular momenta. In Figure0the cor- 
relation of specific angular momenta of different protostellar 
cores with respect to their orientations is shown as a function 
of distance. We calculate the scalar product j ■ j and average 
over cores with similar distances. 

Figure [8^ shows the spatial configuration in model M6k2 
after 15 % of the available material has been accreted onto 
the protostellar cores. These cores form in small aggregates 
with diameters below 0.07 pc. The corresponding Figure 
shows a spatial correlation of the specific angular momenta 
for small distances. The correlation length is approximately 
0.05 pc. Thus, correlation length and cloud size are closely con- 
nected. This can be understood because within one molecular 
cloud clump neighboring cores accrete from the same reser- 
voir of gas and consequently gain similar specific angular mo- 
mentum. In the early phase of accretion we therefore expect 
disks and protostellar outflows of neighboring protostars to be 
closely aUgned. 

Indeed, several examples of parallel disks and out- 
flows have been reported in low-mass, isolated Bok glob- 
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Fig. 6. Average specific angular momentum of protostellar objects in different turbulent environments as a function of the 
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mass bins were considered. The error bars show the standard deviation of < j >. For more clarity the symbols are distributed 
around the corresponding Mach number M (^ —5%) 



ules by ^Froebrich & Scholz ('2003"), Ka mazaki et al ] (|2003!), 
iNisini et aL (.2001.) . and .Saito et al.. Cl995ii) . Alternative expla- 
nations for the alignment of the symmetry axes of young stars 
include density gradients in the pre stellar phase or the presence 
of strong magnetic fields. However. lDuchene & Menardlll2003l) 
found that the disks of T-Tauri stars driving jets or outflows are 
perpendicular to the magnetic field but disks of T-Tauri stars 
without jet are parallel to the field lines. This is very puzzling, 
showing the complexity of the situation that will naturally arise 
in strongly turbulent flows. 

During subsequent accretion the correlation length de- 
creases to values below 0.015 pc (Fig. Eb). This means that 
only close systems maintain correlated (see also Fig.|8j3). This 



has three reasons. First, small systems of embedded cores are 
likely to dissolve quickl y as cl ose encounters lead to ejection 
(e.g., Reipurth & Clarke l200lh . Only close bin aries are able 
to survive for a long time (e.g., Kroupa 1995a b). The corre- 
lation length therefore decreases with time. Second, the same 
turbulent flow that generated a collapsing high-density clump 
in the first place may also disrupt it again before it is fully ac- 
creted. If the clump contains several protostars they will dis- 
perse, again decreasing the correlation. Third, the opposite may 
happen. Turbulence may bring in fresh gas. The protostars are 
then able to continue accretion, but the specific angular mo- 
mentum of the new matter is likely to be quite different from 
the original material. As protostars accrete at different rate, we 
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Fig. 7. The correlation of specific angular momenta of different protostellar cores in model M6k2 with respect to their orientations 
as a function of distance between the cores. As a measure for the correlation the scalar product of different cores was taken and 
averaged over cores that exhibit similar distances between each other High positive values denote co-aligned and high negative 
values denote anti-aligned angular momenta. The three graphs (a)-(c) show three different times at which 15 % (1), 30 % (2) and 
45 % (3) of the available material was accreted. The error bars show the standard deviation of the averaged correlations. 



expect a spread in j to build up, the alignment will disappear. 
Altogether, at later stages of the evolution, we expect the cor- 
relation between the specific angular momenta of close proto- 
stellar objects is disappeared almost completely. This is evident 
in Figure^t- Furthermore, FigurelSj; demonstrates that most of 
the initial subclustering has disappeared by then. 

6. Summary and conclusions 

We studied the rotational properties and time evolution of 
the specific angular momentum of prestellar and protostellar 
cores formed from gravoturbulent fragmentation in numerical 
models of supersonically turbulent, self-gravitating molecular 
clouds. We considered rms Mach numbers ranging from 2 to 
10, and turbulence that is driven on small, intermediate, and 
large scales, as well as one model of collapse from Gaussian 
density fluctuations without any turbulence. Our sample thus 
covers a wide range of properties observed in Galactic star- 
forming regions, however, our main focus lies in typical low- 
to intermediate-mass star-forming regions like p-Ophiuchus or 
Taurus. 

With the appropriate physical scaling, we find the spe- 
cific angular momentum j of prestellar cores in our models, 
i.e. cloud cores as yet without central protostar, to be on av- 
erage {j) - 1 X 10^"cm^s"'. This agrees remarkably w ell 
with observations of cloud cores bv ICaselh etaP ll2002h or 
iGoodman et al. ( 1993 ). Some prestellar cores go into collapse 
to build up stars and stellar systems. The resulting protostel- 
lar objects have on average {j} = 8 x 10'^ cm^ s"' . This is one 
order of magnitude less, and falls into the range observed in G- 



dwarf binaries ("Duauenn ov & Mavorl 19911) . Collapse induced 
by gravoturbulent fragmentation is accompanied by a substan- 
tial loss of specific angular momentum. This is mostly due to 
gravitational torques exerted by the ambient turbulent flow and 
to close encounters occurring when the protostars are embed- 
ded in dense clusters. This eases the "angular momentum prob- 
lem" in star formation without invoking the presence of strong 
magnetic fields. 

The time evolution of j is intimately connected to the 
mass accretion history of a protostellar core. As interstellar 
turbulence and mutual interaction in dense clusters are highly 
stochastic processes, the mass growth of individual protostars 
is unpredictable and can be very complex. In addition, a col- 
lapsing cloud core can fragment further into a binary or higher- 
order multiple or evolve into a protostar with a stable accretion 
disk. It is the ratio of rotational to gravitational energy that 
determines which route the object will take. This is seen in the 
turbulent cloud cores studied here as well as in simulatio ns of 
isolated cores where magnetic fields are important (e.g.. lBossl 
1999). The y6-distribution resulting from gravoturbulent cloud 
fragmentation rep orted here agrees well with y6-values derived 
from observations ( iGoodman et alJl993l) . The average value is 
/3 ~ 0.05. It is important to note, that we find that the distribu- 
tion of /3 stays essentially the same during collapse and accre- 
tion ( see also iBurkert & BodenheimejEoOflt IGoodman et alJ 
Il993h . 

Although the accretion history and thus the evolution of 
the specific angular momentum of a single protostellar object 
is complex, we find a clear correlation between j and mass 
M. This can be interpreted conveniently assuming collapse of 
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Fig. 8. Orientation of the angular momenta {arrows) and spatial distribution of the protostellar cores {diamonds) that formed in 
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an initially uniform density sphere in solid body rotation. Our 
models of gravoturbulent cloud fragmentation are best repre- 
sented by the relation j cc M^^^ . 

When prestellar cores form by compression as part of su- 
personically turbulent flows and then go into collapse and pos- 
sibly break apart into several fragments due to the continu- 
ing perturbation by their turbulent environment, then we ex- 
pect neighboring protostars to have similarly oriented angu- 
lar momentum, at least during their early phases of accretion. 
Star clusters form hierarchically structured, with several young 
stellar objects being embedded in the same clump of molecu- 
lar cloud material. These protostars accrete from one common 
reservoir of gas and consequently gain similar specific angu- 
lar momentum. Their disks and protostellar outflows therefore 
will closely align. Indeed, there are several examples of paral- 
lel disks and outflows seen in low-mass, isolated Bok globules 
/Froebrich & Scholz"2003'; 'Kam azaki et aDl2003t iNisini et alJ 
boo L; 1 aito et al. 1995). During later phases of cluster for- 
mation, the initial substructure becomes erased by dynami- 
cal effects and the correlation between the angular momenta 
of neighboring protostars vanishes. This is in agreement with 



our numerical calculations of gravoturbulent cloud fragmenta- 
tion. They show small groups of close protostellar objects that 
have almost aligned specific angular momenta. As expected, 
the alignment occurs during the early phase of accretion as 
neighboring protostars accrete material from the same region 
with similar angular momentum. During the subsequent evo- 
lution the correlation length decreases. This is either because 
protostellar aggregates disperse, or because infalling new mate- 
rial with different angular momentum becomes distributed un- 
evenly among the protostars. 

Altogether, the process of gravoturbulent fragmentation, 
i.e. the interplay between supersonic turbulence and self- 
gravity of the interstellar gas, constitute an attractive base for a 
unified theory of star formation that is able to explain and re- 
produce many of the observed featur es in Galactic star forming 
regions ( iMac Low & Klesseijl2004h . Our current study con- 
tributes with a detailed analysis of the angular momentum evo- 
lution during collapse. 
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